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Even parity excitations of the nucleon in lattice QCD 
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We study the spectrum of the even parity excitations of the nucleon in quenched lattice QCD. We 
extend our earlier analysis by including an expanded basis of nucleon interpolating fields, increasing 
the physical size of the lattice, including more configurations to enhance statistics and probing closer 
to the chiral limit. With a review of world lattice data, we conclude that there is little evidence of 
the Roper resonance in quenched lattice QCD. 
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I. INTRODUCTION 

One of the long-standing puzzles in baryon spec- 
troscopy has been the low mass of the first positive parity 
excitation of the nucleon, the — Roper resonance, 
or iV*(1440). In constituent (or valence) quark mod- 
els with harmonic oscillator quark-quark potentials, the 
lowest-lying odd parity {J^ = ^ ) state naturally occurs 
below the positive parity radial excitation (with principal 
quantum number iV = 2), whereas in nature the Roper 

resonance is almost 100 MeV below the \~ Af*(1535) 
state. Without fine tuning of parameters, valence quark 
models tend to leave the mass of the Roper resonance too 
high. 

Over the years various suggestions have been made to 
explain this anomaly, ranging from speculations that the 
Roper resonance may be a hybrid baryon state with ex- 
cited glue or a meson-baryon system [Sj], or in terms 
of "breathing modes" of the ground state nucleon Q . To 
understand the nature of the Roper resonance in the con- 
text of QCD, a number of studies have been performed 
recently within lattice QGD. 

The study of excited baryons on the lattice has had a 
relatively short history, although recently there has been 
growing interest in identifying new techniques to isolate 
excited baryons, motivated partly by the experimental 
N* program at Jefferson Lab. The first detailed analysis 
of the positive parity excitation of the nucleon was per- 
formed by Leinweber [5i] using Wilson fermions and an 
operator product expansion spectral ansatz. 

In previous work by the CSSM Lattice Collaboration 
an analysis of the spectrum of octet baryons was per- 
formed using the FLIC fermion action. In each channel, 
a 2 X 2 correlation matrix was used to extract the low- 
lying states. This approach was found to be successful 
in extracting the first excited state of the negative parity 
TV* (1535) state, and in the analysis of the A interpolating 
field. However, the identification of the Roper resonance 
with this correlation matrix remained elusive. 

In the present study we extend the earlier work in sev- 



eral directions, while focusing on the even parity nucleon 
spectrum. In addition, we work with a larger lattice vol- 
ume (2.5 fm compared to 1.95 fm in Refs. 6, 7]), reducing 
finite volume effects and enhancing the statistics. Most 
importantly, we also use an expanded basis of interpola- 
tors compared with that in Ref. Q , with the addition of 
the spin-1/2 projected nucleon interpolator used in the 
calculation of the spin-3/2 hadron mass spectrum 

In the even parity spin-1/2 nucleon channel it is well 
known that the two standard interpolating fields, which 
we label xi S'Hd X2) individually access the ground state 
and an excited state, respectively. The application of a 
2x2 correlation matrix with these interpolators finds 
no evidence of a state with a mass different from those 
that can already be extracted with the two interpola- 
tors individually Furthermore, the extracted excited 
state is found to be too massive to be identified with 
the Roper resonance, and is therefore more likely to have 
stronger overlap with the next even parity excited state 
of the nucleon with mass 1710 MeV — which we denote 
by iV'(1710) (in general we label even parity nucleon ex- 
citations on the lattice by a superscript " ' " ^ and odd 
parity excitations on the lattice by a superscript " * "). 
These findings are consistent with a similar correlation 
matrix analysis by Sasaki et al. @]. 

At the larger quark masses typically used in lattice 
calculations of the spectrum, we expect that the three 
lowest-lying spin-1/2 even parity states are the ground 
state nucleon, the Roper, and the second even parity ex- 
cited state (the A^'(1710)), the latter which appears to 
have strong coupling to the X2 interpolator. One would 
therefore naively expect that the addition of a third nu- 
cleon interpolator to our basis should allow the mass of 
the Roper to be extracted (in quenched lattice QCD). In 
Ref. 7] the spectra of the nucleon and A were analysed, 
including both spin-1/2 and spin-3/2 excited states. For 
the nucleon spectrum a mixed spin-1/2, spin-3/2 interpo- 
lating field (labeled xs) was used. The spin-1/2 projected 
X3 interpolator was found to have good overlap with the 
ground state, and in the present work we use this inter- 
polator as the third interpolating field. 
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In a similar analysis, Brommel et al. used the xi 
and X2 interpolators with the time-component of the xa 
interpolator as a basis for a 3 x 3 correlation matrix anal- 
ysis. Even with the larger basis, Brommel et al. do not 
identify the Roper on the lattice. The difference between 
that study Q and our present analysis is that we con- 
sider the spatial components of the X3 interpolator, with 
spin-1/2 projection as in Ref. 

In Sec. In] we review existing lattice calculations of the 
positive parity excited nucleon spectrum and attempts 
to identify the Roper resonance on the lattice. Our lat- 
tice techniques are outlined in Sec. IIIIl where we firstly 
summarise our simulation parameters and interpolating 
fields. This is followed by a discussion of how to identify 
the spinor indices in which the odd and even parity con- 
tributions to the correlation functions propagate. Our re- 
sults are reported in Sec. lIVl and conclusions summarised 
in Sec. Ivl 



II. EXISTING LATTICE RESULTS 

In this section we review the findings of earlier lattice 
studies of the spin-1/2, even parity nucleon mass spec- 
trum. Figure [T] shows a compilation of recent calculations 
of the mass spectrum in quenched lattice QCD. Because 
the masses of the excited states at small lattice volumes 
are expected to suffer from significant finite volume ef- 
fects, we focus only on those results obtained on lattices 
with a physical size > 2.0 fm. 

In Fig. [T] the studies using Bayesian techniques — 
namely, Sasaki et al. [ll| (pluses) and Mathur et al. 
[13] (crosses) — identify an excited state which is inter- 
preted as the Roper resonance. The physical size of the 
lattice in these analyses is 3.0 fm and 3.2 fm, respec- 
tively, and both use point sources. In both studies the 
mass of the odd parity excited state is found to be consis- 
tent with that of the empirical iV*(1535) resonance. At 
large quark masses, the level ordering of the even and odd 
parity excited states is reversed compared to the physical 
level ordering. However, at small quark masses (ttItt as 
small as 180 MeV) Mathur et al. [12] find that the cor- 
rect empirical ordering appears to be restored. However 
it is suggested by Mathur et al. that greater statistics are 
required to determine if this is true. 

The masses in Fig. [T] which are extracted using a cor- 
relation matrix include the previous study by the CSSM 
Lattice Collaboration [^ (open circles), Brommel et al. 
[1] (squares), Burch et al. flS] (triangles), and Basak 
et al. [13] (diamonds). The mass of the excited state 
extracted with a 3 x 3 correlation matrix in the present 
analysis is also shown (filled circles). 

The previous CSSM work used the xi and X2 in- 
terpolators as a basis in the correlation, with a physical 
lattice size of 1.92 fm. At the fermion source 20 sweeps 
of gauge invariant Gaussian smearing were used, with a 
smearing fraction of 0.7. The findings strongly suggest 
that the mass of the extracted excited state is too large 



to be identified with the Roper resonance. These results 
sit somewhat high because of the relatively small volume 
employed in the analysis. 

A larger basis of operators is considered by Burch et al. 
in Ref. [13|]. That study considers xij X2 and an inter- 
polator equivalent to the temporal component of X3 (up 
to an overall factor of 70). To expand the operator basis, 
two different fermion source and sink Jacobi smearing 
prescriptions, labeled "wide" and "narrow", are consid- 
ered for each quark. The narrow sources have 18 sweeps 
of smearing with a smearing fraction of k = 0.210, while 
the wide sources have 41 sweeps with k — 0.191. The 
physical size of the lattice is 2.38 fm. In the analysis 
of Ref. [1^ the basis of operators is restricted to xi and 
the X3-hke interpolator with three different smearing pre- 
scriptions for each interpolator, making a total of six dif- 
ferent operators. 

In Refs. [13] and 14] Burch et al. argue that in the 



limit of large quark mass, the state corresponding to 
the N' can be identified by comparison with the mass 
extracted with the X2 interpolator. They consequently 
identify the largest mass state with the N' and conclude 
that the lower energy state is therefore the Roper. We 
note that at the larger quark masses the lower energy 
state, identified as the Roper with the Bayesian tech- 
niques in Ref. [12], is similar to the mass of the same 
state identified using a correlation matrix analysis in 
Ref. [l3]- However, the two techniques disagree at the 
smaller quark masses. 

Basak et al. [l3] take advantage of the discrete symme- 
tries on the lattice and identify a large basis of local and 
non-local operators [13] ■ The complete set of local and 
singly displaced non-local sources are used as a basis of 
their correlation matrix, and two distinct excited states 
are identified. The interpretation of Basak et al. [13] is 
that the lower energy excited state corresponds to the 
N'. This state is consistent with the higher energy state 
found by Burch et al. [l3], and our new results presented 
herein. 

The best fit to the excited state data extracted on a 
lattice with physical size ^ 2.5 fm using the correlation 
matrix technique is shown by the solid line in Fig. [T] 
We fit the data for the largest energy state extracted 
by Burch et al. [l3], along with the data from Basak 
et al. [13] and Brommel et al. [3] , and the results of the 
present study. Our results are consistent with the masses 
extracted by Brommel et al. and by Burch et al. at their 
largest quark masses. At the smaller quark masses, the 
highest mass states obtained by Burch et al. and Basak 
et al. lie on either side of the line of best fit. Since these 
two analyses use different gauge and quark actions, it is 
not possible to determine whether the small discrepancy 
between their masses is of any statistical significance. It 
is also likely that the operators considered in each study 
have different couplings to the TV -I- tt in P-wave (and 
equivalently the P-wave N + ij' in quenched QCD) scat- 
tering state. We know that as smaller quark masses are 
approached, the level ordering between the lowest energy 
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(A^O) N 

FIG. 1: Compilation of current lattice calculations of the spectrum of the even parity spin-1/2 excited states of the nucleon. In 
each study the ground state and excited state masses are shown. Burch et al. (closed triangles) report two excited states, one 
as the A'^' the other as a Roper candidate. Basak et al. (diamonds) report the masses of three excited states. Note that the two 
data points with degenerate masses from Basak et al. [l(| have been displaced horizontally for clarity. To aid in understanding 
the results, lines of best fit to the A*'' and ground state masses are included. The solid line is a line of best fit to the N' excited 
state masses extracted in this study, and the studies of Basak et al. , Burch et al. and Brommel et al. . The dashed line is a 
line of best fit to the ground state mass for all of the data shown. 



4 



multi-hadron state and the N' state is reversed on the 
lattice. 



III. LATTICE TECHNIQUES 

The present analysis is based on an ensemble of 396 
gauge-field configurations on a 20^^ x 40 lattice, using 
the mean-field C'(a^)-improved Luscher-Weisz plaquette 
plus rectangle gauge action [l^. The lattice spacing is 
O.f 28 fm, set with the Sommer scale ro = 0.49 fm. For the 
fermion propagators we use the FLIC fermion action [13] , 
which is an C'(a)-improved action with excellent scaling 
properties, providing near continuum results at finite lat- 
tice spacing [18]. A fixed boundary condition in the time 
direction is implemented by setting Ut{x, Nt) = V a; in 
the hopping terms of the fermion action. Periodic bound- 
ary conditions are imposed in the spatial directions. We 
find that the fixed boundary effects are only significant 
after time slice 30 1,19], which is the limit of our analysis 
of the correlation functions presented below. 

We apply 36 sweeps of gauge-invariant Gaussian 
smearing pOj . with smearing fraction a — 0.7, in the 
spatial dimensions at the fermion source (t — 8). Eight 
quark masses are considered in this calculation, providing 
am^ = { 0.540, 0.500, 0.453, 0.400, 0.345, 0.300, 0.242, 
0.197 }. The error analysis is performed by a second- 
order, single-elimination jackknife, with the per de- 
gree of freedom obtained via covariance matrix fits. Fur- 
ther details of the fermion action and simulation param- 
eters can be found in Refs. El] and respec- 
tively. We apply the variational method as discussed in 
Refs. @,|2l|. 

A. Interpolating Fields 

In the previous CSSM study [6| the xi and X2 inter- 
polators were found not to have significant overlap with 



each other. Following the approach by Brommel et al. 
0, we extend this analysis by including the xs inter- 
polator that was used by Zanotti et al. [3] to extract 
spin-3/2 nucleon excited states. Our basis of nucleon in- 
terpolating fields is then given by: 

Xiix) = e"''=(«^»(x)C75rf''(a;))w^(x) , 

X^(x) = e'''^{u^''{x)C^,j'^d\x)h,u%x) . (1) 

In all of our phenomenology we use the Dirac represen- 
tation of the 7-matrices. 



B. Excited baryons on the lattice 

We begin our discussion with a review of how the 
masses of even and odd parity states are extracted from 
the correlation function using the spin-1/2 xi and X2 
interpolators. On the baryon level, the two-point corre- 
lation function in momentum space is: 

git,p) = e-'''-'{0\Txix) m 10) , (2) 

X 

where the interpolator xix) annihilates (creates) baryon 
states to (from) the vacuum. Inserting a complete 
set of intermediate momentum, energy and spin states 
\B,p',s), 

1= \B,p',s){B,p',s\ , (3) 

B,p ',s 

we obtain 



J 



Sit^P) = E E^"*'"^ ^(^) \B,p',s){B,p',s\ m ]0) , (4) 

where the state B has mass Mb and energy Eb ~ \/ M'^ + . The sum over all possible states B with a given 
set of quantum numbers includes a tower of resonances and multi-hadron states created by our interpolators. Using 
x{x) = e*^'^x(0)6~'^'^, where P is the four-momentum operator, we can write: 

g(t,p) = J2 E^^^'"" <0| e^''-"x(0)e-'^-^ \B,p',s){B,p',s\ x(0) ]0) 

s.p ' ,B X 

= E e-^^-'Y.^-"<P-P'^ (0] x(0) \B,p\s){B,p',s\ x(0) ]0) 
= Y e-^^-*%p' (0] x(0) \B,p',s){B,p\s\ x(0) ]0) 

s,p ',B 

^ Ee-^"*E ^(0) \B,p,s){B,p,s\ x(0) ]0) , (5) 

B s 
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where on the last line we make the replacement it —f t for Euclidean time. 

Next we evaluate the matrix elements in Eq. labeling the even and odd parity contributions to the correlation 
function by "+" and "— " , respectively. The overlap of x and x with even and odd parity baryons, such as the nucleon 
for example, can be expressed as: 



X(0) |A^i/2+(p,s)) = Ajv, 



M 



N 



1/2 + 



\ E 



N 



1/2 + 



(0| X(0) |iVi/2-(p,5)) = Xn^ 

{Ni/2+ip,s)\ x(0) |0) = Ajv, 



N 



1/2- 



\ E 



N 



1/2- 



1/2 + 



\ E 



N 



1/2 + 



{Ni/2-iP,s)\ m |0) = -\n,^^ 



I Mm 



E 



N. 



(6) 



1/2- 



where u{p,s) is a Dirac spinor and A(A) are couplings of the interpolators at the sink (source). Note that the four- 
momentum pm , is on-shel] 

^-"1/2 + 

not equal to the adjoint of A. 



momentum Pn-^^^^ is on-shell, with pq = Jjp + Mj^ ^ ^ . Because the fermion source is smeared, the coupling A is 



r 



Substituting the appropriate terms and using the iden- correlation functions, 
tity 



^u(p, s)u(p, s) 



il-p + M) 
2M 



(7) 



the contributions of the even and odd parity terms to the 
correlation function can be written as: 



-E„+t il-PB+ +Mb+) 



Git^P) = VAb+Ab+c 

^ 2£;s+ 



G±(t,0) = tr,p [r±e(i,0) 

= A_B± Xb± exp {-Eb± t) 

^ Ao±Ao± exp(-Afo±t) , 



(9) 



2Ef 



(8) 



where r± = is the parity projection operator at zero 
The masses of states with definite parity can then be momentum [2j|, and the subscripts 0^ label the lowest 
obtained from the spinor trace of the parity projected energy state with the projected quantum numbers. 

I 

For the X3 interpolator the analogue of Eq. (O is given by: 

g^'^it,p) = ^e-^-*^ (0| x^(0) \B,p,s){B,p,s\ x^O) |0) , (10) 

B s 

where now X3 overlaps with both spin-1/2 and spin-3/2 states. In this study we project spin-1/2 states using the spin 
projection operator discussed in Ref. |7|], 



Pi?u (p) = 2 7^71/ + ^ (7 • Pl^^Pi^ + Pt^lul ■ p) 



(11) 



At zero momentum p = (M, 0, 0, 0), relevant to the mass determination, the spin-projection operator has no hypercubic 
lattice artifacts. 

We proceed in our analysis of the X3 interpolator by evaluating the analogue of the matrix element in Eqs. ©. 
Following Zanotti et al. Q, the coefhcient of the spinor is taken to be a linear combination of four- vectors. The require- 
ment that the matrix elements transform as pseudovectors under parity restricts the coefhcients to be proportional 
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to either the four-momentum or the matrix 7''. The analogue of Eqs. ([6]) can therefore be written as: 



^^1/2+ 

IMn , _ 

ix^ (0) |A^i/2- (P, s)) = {aN,,,_ p%^^^_ r)J Hpn,,,- , s), (13) 

V ^1/2- 



(^i/2+(p,s)lx^(0)lO) = ~W-^^ zi(p^,,,+ ,5)75(aiv,,,+ P^^^^^^+^iv,,,+ 7^), (14) 

V -"1/2+ 



(A^i/2-(p,s)|x^(0)|0) = W-^-i^ u(p^^^^_,s)(a^^^^_p^^^^^_+/?Ar^^^_7^), (15) 

where the factors as and Pb denote the couphng strengths of the interpolating field Xs to the baryon B, and similarly 
for the adjoint. Combining these expressions with their respective adjoints, and using Eq. ([7]) for the energy projector, 
the contribution to the correlation function from spin-1/2 states extracted with X3 is: 

B+ "1/2+ 
..^ , 7 • Pm + Mm 

+ T.^-^^-' (^^u.-PK.- ^P^.n-^") ^(«A^./2-^5V,,,_ . (16) 

B- 1/2- 

In this analysis we average over the contributions to the spin-projected correlation function with /i ~ 1-3 and v — 1-3. 
We need not evaluate the (/x, v) = (fc, 0) or (0, k) terms, where k = 1-3, as these do not contribute to the correlation 
function after spin projection at p = 0. 

Finally, we evaluate the two-point function at the hadronic level for the cross correlators (OlTxaXilO) and 
(0|TxiX3|0)j with i = 1,2. It is important to note that as these correlation functions are not Lorentz scalars they 
remain dependent on the representation of the 7-matrices. The respective correlation functions are given by: 

S^dt,P) = E exp(-*p • x) (0 IT xi;{x) x.(0) | 0} , (17) 
Q^^P) = E exp(-*p ■ x) (0 IT x^{x) X^(0) I 0) . (18) 

X 

As for the diagonal correlators discussed above, we proceed by inserting a complete set of states and evaluating the 
resulting matrix elements. For the function 5^'^, we can use Eqs. ([6]), ([7]) and l|12p to write the matrix elements as: 



E^0|^3(0)|A^i/2+(p,s))(A^i/2+(p,s)|x.(0)|0) 



(7 • Pn , i + Mn , ) 

= V/2. i--.^..Pk,.. +/3-V2.^^) ^ ■ (19) 

"1/2 + 

At p = 0, the positive parity states for = 1 propagate in the real part of the (1,2) and (2,1) spinor elements 
of the correlation function. For fi = 2, the positive parity states propagate in the imaginary part of the (1,2) and 
(2, 1) elements, with a relative minus sign. For /i = 3 the positive parity states propagate in the real part of the (1, 1) 
and (2,2) elements, with a relative minus sign, and for ^ = they propagate in the real part of the (1,3) and (2,4) 
elements. 

Similarly, the odd parity contribution to the correlation function is: 



E(0|x^(0)|A^i/2-(p, s)}{N,/,- (p, s)lx.(O)lO) 

S 

(7 • Pn , — ) 

A.,,,_ iaN^^^-P''. + (^N,,.-rh.—^^ ^ . (20) 



N. 



1/2- 



Combining the even and odd parity contributions, we obtain for the "3z" correlation function: 

G!^.(t,P) - E^""-^*^-v2. («^v2.P^v.. +PN,^2.r) ^ /^-^"v^^ +"^"^2.) 

B+ ^1/2+ 
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2£; 



B- "1/2- 

Using the appropriate terms in Eqs. ^ and (I12p in Eq. (fT5)) . the "z3" correlation function can be written: 



1/2 + 



B+ 



2E 



^B-t AjYi/,- 75 



(7 -^7^2- 



(21) 



(22) 



These functions can then be used to relate the appropriate elements of the correlation function to a particular parity. 
To improve our statistics we will average the correlation fun ctions over the spatial components of x-^- 



IV. 



RESULTS 



We begin our analysis of the spectrum by considering 
the 2x2 correlation matrices with xi and Xg, and with 
X2 and Xg. The 2x2 correlation matrix with xi and 
X2 has previously been explored in Ref. . To choose a 
time slice at which to invert the correlation matrix, we 
determine the earliest plateau available to the individual 
interpolators. Our initial time is taken to be one time 
slice earlier than this. If this analysis fails, we invert the 
correlation matrix at one time slice earlier. This algo- 
rithm is discussed in more detail in Ref. [21)] . Throughout 
our correlation matrix analysis we only consider a shift 
of one time slice from this inversion time. 

Proceeding with the 2x2 correlation matrices, the 
eigenvectors for the projection of the xi, Xa correlation 
matrix are obtained from an analysis at time slice t = 14, 
which is two steps back from the onset of the plateau of 
the effective mass extracted with X3 , and six time slices 



after the source. For the X2, X3 correlation matrix, the 
onset of the plateau in the effective mass extracted with 
the X2 interpolator is at t = 11. Eigenvectors for the 
projection of the X2, X3 correlation matrix are therefore 
obtained from an analysis at t = 10. 

The masses of the ground and excited states for each 
correlation matrix are shown in Fig. [21 along with the 
masses extracted with the Xi and X2 interpolators indi- 
vidually. In each case we find that the mass of the ground 
state extracted with the correlation matrix analysis is in 
excellent agreement with the mass of the state extracted 
with the Xi interpolator. Furthermore, the mass of the 
excited state extracted with the correlation matrix agrees 
well with the mass extracted with the X2 interpolator. 
The previous CSSM study 6] showed that the xi inter- 
polator is largely orthogonal to X2- The present calcu- 
lation shows, therefore, that the Xa interpolator has a 
significant overlap with the states accessed by both xi 
and X2- 



With this information in mind we need to determine if the spin-1/2 projected X3 interpolator is a simple linear 
combination of xi and X2- Using the Fierz identity: 



(23) 



where F is one of the matrices in the Dirac representation {1, 75, 7,^, 7,^75, cr^i/|;j>iy} , one can show that: 

Xl^ = ie'^''=(M°'^C7'^u^)7<,7^d'= - le'''"'{u''^Ca''Pu'')a„p-f^'d'' . 
4 8 

Acting on X3 with the spin-1/2 projector from Ref. we can write: 

P^i^Xs = ^7m 71^X3+ ^{l-Plf^Pi^+l^Pt^l-p)X3 
Expanding the combination 71^X3 using Eq. (|24p . we obtain: 

_ !-* ^ ahc f aT /"I a h \ J> jc ^ ahc f aT /"i an b\ J> jc 

7>^X3 = (" C7 M ) 7,y7<T7 " - g£ (" e'er ^'u ) -ivCapl d 



1 



(24) 



(25) 



1 



:75XSR 



(26) 
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FIG. 2: Masses extracted with a 2 x 2 correlation matrix with xi a-nd Xa (left), and xa and Xa (right). For comparison the 
masses extracted with the xi and X2 interpolators are also shown. The data correspond to m,r — 830 (rightmost points), 770, 
700, 616, 530 and 460 MeV (leftmost points). 



where we identify XSR from Rcf. 5] as the Unear combination 2{x2 ^ Xi)- Substituting this result in Eq. (|24[) . we find: 

(27) 

At zero momentum the spin-1/2 projected Xs can then be written: 



Pkl. Xs = gTsTfclXi - X2) + g7o7feX3 



(28) 



r 



Since we consider the spatial components of the spin pro- 
jected correlation function, fi ~ 1-3, the only new infor- 
mation in the Xs interpolator is from X31 i-S- the time 
component of the X3 interpolator used by Brommel et al. 
0. Thus our analysis serves as an independent check of 
Ref. 0. 

The Fierz transformation of the X3 interpolator allows 
us to verify that this has a strong overlap with both the 
Xi and X2 fields. The question remains whether there 
is sufficient additional information in the X3 interpolator 
to extract a second excited state. We proceed, therefore, 
with a 3 X 3 correlation matrix analysis with xii X2 a-nd 
X3. Eigenvectors for the projection of the correlation 
matrix are obtained from an analysis at t = 10, one step 
back from the onset of the plateau in the effective mass 
extracted with the X2 interpolator. 

The extracted masses are shown in Fig. [3l along with 
the masses extracted with the xij X2 and X3 interpola- 
tors individually. All three determinations of the ground 
state mass are found to be in excellent agreement, as are 
the masses of the excited state extracted with the correla- 
tion matrix, and the X2 interpolator. At the two smaller 
quark masses we only have sufficent statistics to fit two 
of the three diagonal elements of the projected correla- 



tion matrix. The masses extracted with the xi and X2 
interpolators individually are shown in Table. [D 



TABLE I: The mass of the nucleon extracted with xi and the 
even-parity spin-1/2 state excited state extracted with X2. 
For comparison we show the mass of the even-parity spin-3/2 
state extracted with X3- 





mjv 


mjv' 


3/2+ 


0.541(1) 


1.061(3) 


1.633(14) 


1.655(9) 


0.501(1) 


1.012(4) 


1.609(18) 


1.624(10) 


0.454(1) 


0.955(5) 


1.589(21) 


1.590(11) 


0.401(1) 


0.895(6) 


1.576(25) 


1.557(13) 


0.347(2) 


0.836(9) 


1.581(34) 


1.531(13) 


0.302(2) 


0.790(12) 




1.515(16) 


0.241(4) 


0.755(13) 




1.515(19) 


0.197(5) 


0.745(14) 




1.545(26) 



Our results show no evidence of a Roper-like even par- 
ity excited state, which suggests that the couplings of the 
interpolators used in this study to such a state must be 
either small or zero. On the other hand, using Bayesian 
techniques the analyses in Refs. [ll|i do report a low- 
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FIG. 3: As in Fig.[2j but for the 3x3 correlation matrix with 
the xi, X2 and Xs interpolators. 

lying state. If the Roper resonance exists in quenched 
QCD near its experimental value, then it would appear 
that our smeared sources have unfortunately poor over- 
lap with this state. However we note that Brommel et al. 
P|, Burch et al. [l^, and Basak et al. [l3| all report re- 
sults similar to ours with different interpolating fields. 
It is difficult to understand, therefore, if this state does 
actually exist, why it would not be seen in any of these 
analyses. 

In Fig. [4] we enlarge a portion of Fig. [U adding the 
energies of the non-interacting P-wave N + tt for each 
study, and the mass of the even parity spin-3/2 nucleoli 
calculated here and in Ref. [T3|. We show the mass of 
the spin-3/2 state as a guide because it has a mass of 
1720 MeV, similar to the mass of the N'. Further, at 
the larger quark masses it is reasonable to expect that 
the hyperfine splitting between the spin-1/2 and spin- 
3/2 states should be small. Recall that the energy of the 
P-wave N + TT threshold state will be larger than the at- 
tractive scattering state mass due to finite volume effects 
and that a stable state is likely to appear at sufficiently 
heav y qu ark masses [l^, [2l| . Note that for the data from 
Ref. [l3| labeled as spin-3/2, only one of the three degen- 
erate masses corresponds to a spin-3/2 state, with the 
others corresponding to higher spin states. 

At the three largest quark masses, the mass of the 
excited state calculated in this study is consistent with 
the highest-lying excited state calculated by Burch et al. 
and with the spin-3/2 state. We also reproduce the mass 
of the excited state calculated by Brommel et al. . As 
we approach smaller quark masses, the mass of the ob- 
served excited state becomes larger than the energy of 
the non-interacting P-wave N + Tr state. In this case the 
reported masses will include some contamination from 
the Ntt threshold, although the spectral strength may 
be relatively weak. Even with this low-mass contami- 



nation there is no evidence of a state approaching 1440 
MeV. 



V. CONCLUSION 

In this study we have attempted to extract the mass of 
the first even parity excitation of the nucleon in quenched 
lattice QCD, using a basis of three nucleon interpolating 
fields. We have demonstrated that the ground and second 
even parity excited states of the nucleon can be extracted 
from a 2 X 2 correlation matrix, with either the xi and 
X3, or X2 and X3 interpolating fields. With the use of 
the Fierz identity, we showed that the spin-1/2 projected 
X3 interpolator does indeed have a large overlap with the 
Xi and X2 interpolators, and with Xs- 

Extending the analysis to a 3 x 3 correlation matrix, 
we found no evidence that our interpolators overlap with 
the Roper resonance. Our results are in accord with other 
correlation matrix based analyses using similar interpo- 
lating field constructs. The absence of a low-lying ex- 
citation in any of these analyses raises the question of 
how this state is seen in Bayesian analyses. Since the 
multiple-operator correlation matrix approach is an es- 
tablished method for extracting excited state masses, fur- 
ther careful examination of the Bayesian analyses would 
seem appropriate. 

By comparison with the mass of the even parity spin- 
3/2 state at large quark masses we identify the excited 
state extracted in our correlation matrix as the A^'(1710). 
However we recognize that approaching the chiral regime 
the level ordering of the N' and the P-wave N + n (and 
P-wave N + rj' in quenched QCD) is reversed. This ambi- 
guity makes an analysis to discriminate scattering states 
central to future calculations of the spectrum at light 
quark masses. 

On a final note we emphasize the need to bring 
these advanced analysis techniques to bare on dynami- 
cal fermion configurations. There is now ample evidence 
that this is likely to be the key missing ingredient in cre- 
ating and observing the Roper resonance in lattice QCD 
simulations. 
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